High-fidelity quantum operations on superconducting qubits in the presence of noise 
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We present a scheme for implementing quantum operations with superconducting qubits. Our 
approach uses a "coupler" qubit to mediate a controllable, secular interaction between "data" qubits, 
pulse sequences which strongly mitigate the effects of 1/f flux noise, and a high-Q resonator-based 
local memory. We develop a Monte-Carlo simulation technique capable of describing arbitrary noise- 
induced dephasing and decay, and demonstrate in this system a set of universal gate operations with 
0(10"^) error probabilities in the presence of experimentally measured levels of 1/f noise. We then 
add relaxation and quantify the decay times required to maintain this error level. 

PACS numbers: 02.70.Ss, 03.65.Yz, 03.67.Ac, 03.67.Lx, 03.67.Pp, 85.25.Am, 85.25.Cp 



Superconducting qubits are promising building blocks 
for quantum computers [ij. Realizing this promise 
requires a qubit coupling scheme and associated con- 
trol pulse sequences capable of realizing gate operations 
with error probabilities low enough to achieve fault- 
tolerance Such a control scheme should ideally per- 
mit a practical, time-efficient implementation of universal 
single- and two-qubit operations; incorporate switchable 
coupling between qubits; and be insensitive both to de- 
coherence and fabrication variations of the qubits. 

Although several proposed [1, 0, [E Hi and experimen- 
tal qubit coupling schemes already exist, none yet 
combine all of the above features. In this Letter, we 
present a scheme which aims to do this. Our approach 
utilizes a "coupler" qubit to mediate interaction between 
"data" qubits [J] , pulse sequences incorporating spin- 
echoes to suppress the effects of 1/f flux noise [sJm|> and 
local memory in the form of high-Q resonators 10]. We 
introduce a Monte-Carlo technique capable of simulating 
gate operations in the presence of noise and spontaneous 
decay, which we then use to show that C(10~^) error 
probabilities (comparable to certain predicted thresholds 
for fault-tolerance 0) can be achieved in the presence 
of experimentally measured levels of flux noise. Finally, 
we quantify the qubit excited-state lifetimes and control 
signal requirements necessary to achieve this. 

We focus here on the flux qubit 11], although much 
of what follows is more generally applicable. In most 
work to date, this consists of a single superconduct- 
ing loop interrupted by several Josephson junctions [Fig. 
[IJb)]. In this work, however, we consider the double-loop 
four-junction version [Fig. [TJc)]. When ^1,^2 ~ ^,0 
(^o = h/2e is the flux quantum), the qubit can be de- 
scribed by the approximate two-level Hamiltonian Hqb — 
+ A(T^]. The eigenstates of ct^ are "persistent cur- 
rent" (PC) states with energies ±e/2 = ±/p($i - ■^), 
where currents ±/p circulate around the <I>i loop. The pa- 
rameters e and A are tunable by the fluxes $e = $1 + ^ 
and $A = ^'2 [Figs. Point A in Fig. [T^a), 

where the PC states are mixed to produce an avoided 
crossing with ir^-like energy eigenstates, is known as the 
"degeneracy point". It has the desirable feature that 



FIG. 1: (color online) Flux qubits and secular vs. nonsecu- 
lar coupling, (a) Lowest-two energy levels for a flux qubit. 
Schematic of (b) single-loop 3-JJ and (c) double-loop 4-JJ 
flux qubits. Relative junction areas have ratios a and fi [ll| . 
Single-loop qubits have only one control flux "5^ and are re- 
stricted to the line AC in (a). Double- loop qubits have two 
control fluxes: $e = $1 + and <1>a = "I>2. (d) Energy 
levels for two flux qubits at degeneracy, with nonsecular cou- 
pling. Level shifts result from second-order mixing of two- 
qubit eigenstates (dashed arrows), (e) Energy levels for sec- 
ular coupling. A conditional level shift allows CNOT to be 
implemented with a single microwave transition (solid arrow) , 
and SWAP with a stimulated Raman transition (by driving 
a 2n pulse on each qubit, corresponding to a tt pulse on the 
Raman transition), (f) possible layout of a gradiometric 4-JJ 
qubit, for which "I>i and $2 can be accessed independently. 
The A control-loop orientation induces $1 and $2 in the - 
0.5:1 ratio required to adjust A without affecting e. 
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(where the qubit energy split- 
so that decoherence due to flux 
Existing qubit coupling schemes 



noise is minimized 
which are insensitive to decoherence rely on biasing at 
this point 

There are, however, two disadvantages to working at 
degeneracy: first, E is then fixed by fabrication and may 
vary significantly between qubits, requiring individually 
tuned microwaves; second, the degenerate qubit is first- 
order insensitive to flux from other qubits, so that indue- 
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FIG. 2: (color online) Possible scheme for universal two-qubit 
operations, (a) A 4-JJ coupler qubit mediates the interac- 
tion between two 4-JJ data qubits (Ql, Q2). Coupling is 
switched oflt by tuning ec ^ Ac. (b) An inductive interac- 
tion mixes the coupler levels (dashed arrows), resulting in a 
conditional frequency shift 5vc in the computational subspace 
C2 [see Fig. [TJe)] (c) Increasing the qubit A adiabati- 

cally through the avoided crossing transfers the qubit state to 
the resonator [lo| . 



tively coupling <i>i -loops, as in previous works [1, 13, 
has no effect on the system to first order. This "nonsec- 
ular" coupHng can be written in the form H'^ = J^o-^ctI 
(the subscripts indicate qubits 1 and 2), and its leading- 
order effect at degeneracy is a second-order mixing [Fig. 
[Ijd)]. This mixing can be exploited using double reso- 
nance or parametric driving 0, Q to enable, for ex- 
ample, the two-qubit VSWAP gate, even if the qubits 
have different splittings [3|, U, |5|; however, these schemes 
require different microwave frequencies for each qubit or 
qubit pair. Furthermore, the Controlled-NOT (CNOT) 
gate used ubiquitously in quantum circuits must then 
be constructed from several vSWAP and single-qubit 
gates [T2j . 

Far from degeneracy (e ^ A), however, the energy 
eigenstates are approximately PC states, and com- 
mutes with Hqh, giving a first-order level shift [Fig. [TJe)] 
which hfts the degeneracy between the | ||) | ti) and 
I it) ^ I TT) transitions by the conditional frequency 
shift hSvc — Jc- "^^i^ "secular" coupling allows CNOT 
to be driven with a single microwave pulse (Fig. [IJe), 
solid arrow), and SWAP with a stimulated Raman tran- 
sition (dashed arrows). However, the qubits are now sen- 
sitive to flux noise due to the large slope dE/d^^ = 2/p. 

Based on these observations, we consider the system 
shown in Fig. dja) consisting of two two-loop "data" 
qubits, whose coupling is mediated by a third "coupler" 
qubit in a manner similar to that used in Ref. J, for single- 
loop qubits. Here, however, we inductively couple to <i>A 
of each data qubit rather than $1 (as in previous works), 
with the Hamiltonian: Hq = Jc(Tq(i7^-|-(T2), where Jq ~ 
MicIpdAi/d^s and Mic is the mutual inductance be- 
tween data and coupler qubits [l^. When ei = €2 = ec — 
0, this configuration produces the level structure shown 
in Fig.[2fb) [14]. The lower four of these eight eigenstates 
(hereafter the two-qubit "computational subspace" orC2) 



exhibit an effective two-qubit conditional frequency shift 
hSvc = Ac{l - y^l + (2Jc/Ac)^} equivalent to that 
shown in Fig. [ije) [15]. In addition, this shift can be 
turned off by tuning ec/Ac ^ 1, so that the coupler 
eigenstates become approximately PC states (and eigen- 
states of He)] this suppresses the mixing shown in Fig. 
Mh), and with it Si^c by the factor (Ac/cc)^- The 
speed of this switching is limited only by nonadiabatic 
excitation of the coupler qubit [l6| . Finally, the sensitiv- 
ity of each data qubit to flux noise in this configuration 
is given by dA/d^A — Jc/MIp {1^ is the coupler per- 
sistent current), which can be reduced substantially by 
using larger junctions in the coupler qubit to increase Ip. 

The final ingredient shown in Fig. [UJa) are the 
transmission-line resonators coupled to the $e-loop of 
each qubit so that its state can be adiabatically trans- 
ferred into a resonator state [13] [Fig. He)]. These res- 
onators provide a potentially long-lived memory in which 
qubit states can be stored between manipulations. 

To simulate our system, we integrate the time- 
dependent Schrodinger equation with the Hamiltonian: 



+ [^^-^rV^J- + + a,,)a^] + He (1) 

where i E {1,2,C}, j G {1,2}, and we have added mi- 
crowave fields illuminating the loop of each qubit. 
We treat the qubits as two-level systems and truncate 
the ladder of resonator photon states at n = 1 [l^ . As 
fixed input parameters, we assume resonator frequencies 
fri,r2 = 9,11 GHz (they may be identical in general), 
qubit-resonator coupling Jrcs/h = 0.75 GHz, and qubit- 
qubit coupling Jc/h = 0.3 GHz 

We have simulated Hadamard (H), CNOT, and SWAP 
gates. As an example, the CNOT pulse sequence is 
shown in Fig. ^a,). Both qubits begin with states 
stored in their respective resonators (Ai,2 high), and 
with the coupler "off' (ec/Ac ^1). A downward 
ramp is applied to both A.^ to retrieve the qubits' states 
(only one A^ would be ramped for H). These ramps 
are optimized to maintain a constant non-adiabaticity 
as the avoided crossing in Fig. [2jc) is traversed. Next, 
the qubits are manipulated with single-qubit microwave 
rotations and conditional phase (C<i>) gates produced 
by transiently pulsing the coupler "on" (ec/Ac = 0) 
using linear flux ramps [T^]. The CNOT pulse sequence 
(with bit 1 the target qubit) shown in Fig. ^a,) is: 
i?i,2(7r)i?f(f)C-(f)i?f(f)C-(f)i?i,2(7r), where R;,^{e) 
is a rotation of angle 9 around axis v on the I31och 
sphere for qubits p and q, Cif) = C(f )i?i,2(7r)C(f ), 
and we have suppressed the free evolution inter- 
vals for clarity. The 7r-pulses induce spin-echoes 
at the center of each 7r/2-pulse and at the end of 
the operation. H and SWAP can be constructed 
similarly: H = Ri{n)Rf {^)Ri{n), and SWAP = 
RMRU^)C%^)Rl,{^)C'^{^)Rl2{^)C%^)RM. 
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FIG. 3: (color online) Simulation of gate operations in the presence of noise, (a) CNOT pulse sequence. Data qubits are ramped 
to Ai,2 =5,6 GHz during operations, e.g., point X [Fig. [ija)] (Ai,2 are chosen to be different here only for clarity); the coupler 
has (ec, Ac) = (17, 0.2) GHz and (0, 4) GHz in the decoupled and coupled states, respectively fl^ ]. Pulsed microwaves generate 
TT and 7r/2 rotations; Ctt/A operations are generated by coupling for a controlled duration. Spin-echoes are generated at times 
indicated by dashed lines using vr- pulses, (b) Simulation of a single qubit including 1/f noise and spontaneous decay (see text), 
(c) Decay rates as a function of e/A. Solid symbols show the fitted exponential decay rate Fe of P22 (constant by construction). 
Open symbols show the decay rate Fg, extracted by fitting |pi2| to the approximate expression: exp[— Fcf — (Fgt)^] with Fo 
fixed at the fitted value. The solid red line is based on the theory in Ref. [2l| with no free parameters, (d) maximum gate error 
£ for H (triangles) and CNOT (circles) gates with Ti — ^ 00 and noise spectral density S{lli) = At^/uj, vs. A^i. Open symbols: 
same operations without spin-echo compensation. Vertical dashed line: experimentally measured noise level of 1 fi^o / \/Hz at 
1 Hz 3. (e) Maximum error probabilities at that noise level versus the qubits' excited state lifetime. In both (d) and (e), (+) 
and (x) show £ for CNOT and H, respectively, with noise and decay on the coupler qubit only. 



Finally, the qubit states are transferred back into the 
resonators, and the final relative phase of each qubit is 
adjusted by varying the time at which this occurs. 

All microwave pulses have a fixed amplitude (with 100- 
ps rise/fall times), so that pulse area is controlled solely 
by their duration. This is crucial for producing effec- 
tive spin-echoes with short, intense pulses at a fixed fre- 
quency, due to the Bloch-Siegert shift [l3|. This fre- 
quency shift is amplitude dependent, so the microwave 
amplitude must be kept fixed in order to compensate for 
it with a fixed frequency offset (here, ^--^5 MHz). 

The use of single-qubit rotations and C$ gates instead 
of the frequency-resolved pulses shown in Fig. [IJe) has 
several advantages: first, it requires only a single mi- 
crowave frequency for all operations (assuming all qubits 
are tuned to the same frequency using $a); second, 
it does not require accurately resolving two frequencies 
spaced by Svc using long microwave pulses. Finally, all 
gate parameters can be adjusted solely by varying the 
timing of the microwave and flux pulse edges. These 
characteristics may also be of interest for other types of 
qubits. 

We simulate the effect of decoherence due to flux noise 
using a Monte-Carlo technique. Noise is added to each 
Ai and Ci in eq. [1] by superposing an independent, dis- 
crete set of random two-level "telegraph" signals [l9| . 
which produce a 1/f average spectral density. For each 
Monte-carlo iteration (all of which can be run in paral- 



lel), a new set of noise functions is generated (one for 
each e. A). From these we obtain a sequence of N time 
intervals {Sti} (typically N ^ several hundred) whose 
union is the full gate time, and over each of which all 
fluctuators are constant. We then integrate the 32 cou- 
pled Schrodinger equations sequentially over these inter- 
vals. The final states for all iterations are converted to 
density matrices and then averaged, corresponding to a 
statistical mixture of the different noise realizations. 

Spontaneous decay is added as follows. At the 
end of each Sri = ti — we construct the den- 

sity matrix pi and apply amplitude damping to the 
qubits: pj = S'^{-fi)S'^{-fi)S^{-fi)p^ where S^{-i)p = 
Efct^^(7)pC/^'(7)nand: 

Here, 7^ = 1 — exp[— (^Ti/Tj'] is the probability of decay 
for qubit j during Sri , Tl is the excited-state decay time 
for qubit j, and transforms from the instantaneous 
energy eigenbasis of the isolated qubit j at time ti (de- 
cay is assumed to occur in this basis) to the fixed basis 
of the simulation [20] . To return to state- vector space at 
the end of each Jr^, we first write pi = J2k ^i\''Pi){''Pi\- 
In principle, we could then pick one normalized ji/'f) 
at random according to the probabilities A*^; however, 
since < 0(7^) for A: 1 (with A*^ in decreasing or- 
der for given ?), convergence of the Monte-Carlo sum 
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would require many iterations to sample any of these 
unlikely events. Instead, we use an expansion procedure. 
For each noise realization, we run two simulations; for 
the first, we choose {ipl) at every time ti g {<i • • -iAr}. 
This corresponds to the leading term in a weighted av- 
erage, having probabihty Vq = Hi -^1 ^'(1 — tN/Ti). 
Next, we use the {A^} obtained from simulation of the 
first term to obtain approximate relative probabilities 
■p(rn, n) for the next order terms, sequences of eigen- 
vectors of the form {ipl,--- - ■ ■ ,V']v} 
with m 1 (exactly one unlikely event), as: V{rn,n) « 
\\ - ■ ■ A^_jA™ (we neglect the possibility of two or more 
unlikely events); we then randomly pick one of these 
terms according to normalized probability distribution: 
P(m, n)/ ^jjj ^ 7^(771, n), and simulate it. Finally, we 
combine the density matrices for these two simulations 
with the relative weights Vo and Ylm n ^("^, n) w 1—Vo- 

Figures EKb) and (c) show a test of this method on a 
single qubit initially in the state (| T) + I i))/'s/2, under 
the influence of both 1/f noise and decay. As expected, 
the excited state decays exponentially, while the coher- 
ence does not [Fig. [SUb)]. Decay rates can be extracted 
from these data, and are shown in Fig. \Elc) as a func- 
tion of e/A. The solid line is a prediction, with no free 
parameters, using the theory from Ref. [2ll . 

To characterize a gate operation, we run the simulation 
(in parallel) for each of a set of chosen initial states {(pg) 
spanning C2. From these outputs we can reconstruct the 
superoperator iSop for the gate (acting on density matri- 
ces in C2), which can be represented as a 16 x 16 matrix 
acting on "vectors" obtained by stacking the columns 
of the 4x4 density matrix (2^ . For decoherence only 
(no decay), each iteration is a linear operation, and only 
four |0o) are required. With decay, however, this is no 
longer true due to the projection and rcnormalization. 



and sixteen |(/>q) are required. 

We calculate the maximum error probability: £ = 
1 — min [(5'|5o""^5opp|'I')] , where So is the superopera- 
tor corresponding to the desired operation, and maxi- 
mization is over normalized state vectors |^) in C2. Fig- 
ure [31[d) shows £ for the CNOT (circles) and Hadamard 
(triangles) gates, both with spin-echo compensation [fig. 
[3lja)] (filled symbols) and without it (open symbols), as a 
function of 1/ f- noise amplitude A,^ , for Ti — > 00 [23, [13] ■ 
The spin-echo compensation strongly suppresses errors 
due to 1/f noise, keeping £ at or below 0(10^^) up to 
the noise level observed in Ref. d, indicated by the verti- 
cal dashed line. Figure [31[e) shows £ for this same noise 
level as a function of the qubits' excited state lifetime, 
which would evidently have to be ~ ms for £ ~ 10~^. 
Finally, (-I-) and (x) show £ with noise and decay added 
only to the coupler qubit, illustrating that although the 
extra qubit does in principle expose the system to more 
noise, the present configuration is insensitive to it. 

Although £ can be used to set bounds on fault- 
tolerance thresholds for quantum circuits, the iSop gener- 
ated by these simulations give much more: a full statis- 
tical description of the errors that occur for an arbitrary 
input state. This may provide a means to compute more 
accurate thresholds for real quantum circuits (assuming 
that noise correlations between successive operations can 
be neglected). Furthermore, it may enable the design 
of optimized error-correction schemes and/or gate oper- 
ation protocols which target a particular noise source, 
enabling higher thresholds for fault-tolerance. 

This work is sponsored by the United States Air Force 
under Air Force Contract #FA8721-05-C-0002. Opin- 
ions, interpretations, recommendations and conclusions 
are those of the authors and are not necessarily endorsed 
by the United States Government. 
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